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ABSTRACT 

The continuum spectrum of the planetary nebula Hf 2-2 close to the Balmer discon- 
tinuity is modeled in the context of the long standing problem of the abundance and 
temperature discrepancy found when analyzing optical recombination lines and col- 
lisionally excited forbidden lines in nebulae. Models are constructed using single and 
double Maxwell-Boltzmann distributions as well as K-distributions for the energies of 
the free electrons. New results for the necessary continuum and line emission coeffi- 
cients are presented calculated with ^-distributed energies. The best fit to the observed 
continuum spectrum is found to be a model comprising two components with dramat- 
ically different temperatures and with a Maxwell-Boltzmann distribution of electron 
energies. On the basis of a x 2 analysis, this model is strongly favored over a model 
with ^-distributed electron energies. 
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- nebular physics - optical recombination lines - collisionally excited lines Balmer 
continuum - abundance discrepancy - electron temperature - electron distribution - 
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1 INTRODUCTION 


The long standing problem in nebular physics related to 
the abundance and temperature discrepancy between the 
results obtained from optical recombination lines (ORL) 
and collisionally excited lines (CEL) has been visited by 
many researchers in the last few decades. Several explana- 
tions for the general trend of obtaining higher abundances 
and lower temperatures from ORL than from CEL of the 
same ions have been put forward. One explanation is a 
multi-component model of planetary nebulae in which low 
temperature, high metallicity components (clumps) produce 
the ORL and are embedded in a high temperature, rela- 
tively low metallicity, hydrogen-rich component that pro- 


duces the CEL (Liu et al 2000). It has been suggested re- 


cently (Nicholls et al 2012) that this discrepancy may be 


largely explained by the assumption of a different elec- 
tron energy distribution, a K-distribution, rather than the 
Maxwell-Boltzmann (MB) distribution which is tradition- 
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ally accepted as the dominant distribution in the low density 
plasma found in the planetary nebulae. 

The debate about the electron energy distribution in 
astronomical objects, including planetary nebulae, is rela- 
tively old and dates back to at least the 1940s when Hagi- 
hara H agihara| ( 1944| suggested that the distribution of free 
electrons in gaseous assemblies deviates considerably from 
the MB distribution. This was refuted by Bohrn & Aller 
( 1947]) who, on the basis of a detailed quantitative balance 
analysis, concluded that any deviation from the Maxwellian 
equilibrium distribution is very small. The essence of their 
argument is that for typical planetary nebulae conditions, 
the thermalization process of elastic collisions is by far the 
most frequent event and typically occurs once every second, 
while other processes that shift the system from its ther- 
modynamic equilibrium, like inelastic scattering with other 
ions that leads to metastable excitation or recapture, oc- 
cur at much larger time scales estimated to be months or 
even years. Bohm and Aller also indicated the significance 
of any possible deviation from a Maxwellian distribution on 
derived elemental abundances. However, nobody seems to 
have considered non-MB electron distributions as a possible 


© 0000 RAS 


2 P.J. Storey & Taha Sochi 


solution for the ORL/CEL discrepancy problem until the 
recent proposal of Nicholls et al (20121. In the light of this 
suggestion, it seems appropriate to ask whether there is any 
empirical evidence from observation for departures from MB 
energy distributions in nebulae. 

A preliminary investigation of this possibility has been 
given recently by |Sochi| ( |2012| ) and |Storey fe Sochi| ( |2013| ) 
who proposed a method based on using atomic dielectronic 
recombination (DR) theoretical data in conjunction with 
observational line fluxes to sample the free electron dis- 
tribution and compare to the theoretical distributions, i.e. 
Maxwell-Boltzmann and re. For all the objects studied by 
|Storey fe Sochi| ([2013 ) the limited observational data ana- 
lyzed were found to be better described by a single MB dis- 
tribution than a re-distribution but the uncertainties were 
sufficient that a re-distribution could not be conclusively ex- 
cluded. A two-component model incorporating low temper- 
ature material was also found to not fit the data as well as 
a single MB model. 

The current paper approaches this problem by model- 
ing the magnitude of the Balmer discontinuity and nearby 
Balmer continuum of the planetary nebula Hf 2-2. This neb- 
ula is an extreme example of the ORL/CEL discrepancy 
trend since its abundance discrepancy factor (ADF), defined 
as the ratio of ORL abundance to CEL abundance, reaches 


an exceptionally high value between 68-84 l 

Liu|2003; Wes- 

son et al 2003. Zhang et al - 2004 Liu et al 2004; Wesson 

& Liu||2004 

Liu et a/||2006 Wesson et a/f2008 

). An excep- 


tionally high temperature discrepancy between the CEL and 
ORL results, which differ by about an order of magnitude, 
has also been obtained for this nebula. 

The shape of the Balmer continuum is principally de- 
termined by a convolution of the n = 2 hydrogen recombina- 
tion cross-section and the free electron energy distribution. 
The shape of the continuum therefore offers the possibility 
of gaining information about the free electron energy dis- 
tribution. We shall see below that in the case of Hf 2-2 the 
continuum shape is primarily determined by the free electron 
energy distribution. In section[2]we outline the properties of 
Hf 2-2 and in section [3] we deal with the theory describing 
the continuum processes and new results relating to recom- 
bination using a re-distribution. The fitting procedure and 
the resulting fits to the spectral data with MB and re dis- 
tributions are described in section [4] and our conclusions in 
section [5] 


2 Hf 2-2 


Hf 2-2 is a faint southern highly symmetric low density plan- 
etary nebula with a central cavity and interior disc-shaped 
multi-shell structure. Its spectrum includes carbon, nitro- 
gen, oxygen, neon, sulphur and argon lines as well as helium. 
Hf 2-2 seems to have a high abundance in He and C with 
a high C/O ratio and an exceptionally strong C II A4267 A 
signature. The nebula, which is in the galactic bulge at a 
distance of about 4. 0-4. 5 kpc, has a spectrally-varying close 


binary central system with a period of about 0.40 day ( Bond 


2000 1 . Hf 2-2 has common spectral features with old novae 
like DQ Her ( Cahn &; Kaler|1971 Maciel 1984 Kaler|19~88 


Liu||2003 Liu et a(||2006 Schaub et aZ||2012 l, and therefore 
being a planetary nebula may be disputed. 



Wavelength [Angstrom] 

Figure 1 . The observed Hf 2-2 spectrum (solid) and the con- 
tinuum fit with a single Maxwell-Boltzmann distribution with 
T = 8800 K (dashed). The fit was optimized on the two longer 
wavelength of the three wavelength segments (shown in the figure 
as horizontal lines): 3585-3645, 3780-3790 and 4180-4230 A. 


This unusual nebula is marked by a number of excep- 
tional features. Hf 2-2 exhibits a very high ADF, possibly 
the highest recorded for a planetary nebula, with very strong 
ORL emission. Hf 2-2 also has the very large C/O abun- 
dance ratio of about 19 (Patriarchi & Perinotto 1994). A 
third unusual feature is a very large Balmer jump associ- 
ated with a rapid fall in the Balmer continuum intensity 
towards shorter wavelengths. An estimation based on the 
magnitude of the Balmer discontinuity indicates an electron 
temperature of about 780-1000 K in sharp contrast to the 
forbidden lines estimation of about 8800 K from an [O ill] 


line (Liu 2003 Liu et al 

2004 Zhang et al 2004 Liu et al\ 

2006 Mcnabb et a([2012 

. The deduced Balmer discontinu- 


ity temperature is one of the lowest observed for a planetary 
nebula. These findings led Liu et al ( 2006 ) to conclude that 
a multi-component model of cold hydrogen-deficient knots 
embedded into a metal-poor nebula may be inevitable to 
explain these exceptional observations. 

With regard to ORL electron temperature derivations 
for Hf 2-2, Bastin (2006) obtained a mean electron temper- 


ature of < 600 K from C II lines, while Liu (2006) derived 
an electron temperature of about 630 K and Mcnabb et al\ 


12012) a temperature of about 3160 K from O n lines. As 
for the stellar temperature of this nebula, it is estimated by 
Liu et ai| ( |2006| ) between 50000-67000 K, while |Schaub ~ei| 


al (20121 set the temperature of the two stars (assuming a 
binary system which seems to be largely accepted) in their 
model to 67000 K and 7500 K. 

The observational data of Hf 2-2 which is used in the 
present paper were obtained in 2001 by a Boiler and Cliivens 
long-slit spectrograph mounted on the ground-based 1.52 m 
European Southern Observatory telescope located in La 
Silla Chile. More details about this data set can be found 


in Liu et al (2006). The observed spectrum is shown in Fig- 
ure^ with the measured flux normalized to Balmer Hll, an 
apparently unblended line close to the Balmer discontinuity. 
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3 THEORY 

To model the magnitude of the Balmer discontinuity and 
the shape of the Balmer continuum we consider the contri- 
butions to the continuum spectrum from recombination of 
H + and He + , H 2-photon emission and the underlying scat- 
tered stellar continuum. Recombination of He 2+ is neglected 
due to the very low He 2+ /H + fraction in Hf 2-2 which is es- 
timated to be about 0.002 ( |Zhang et af|2004 Liu et a^|2006 1. 

The recombination of an atomic ion X + with an elec- 
tron, e~ of energy E to a state X* of the recombined ion, 


X+ + e~ -> X* + hv, 


(1) 


gives rise to continuous emission with emission coefficient 
(energy emitted per unit time per unit frequency and per 
unit particle densities) 


7*M = 


h 3 a 3 c w* ( hu\ 3 ( R \ 2 


327r 2 m e aQ a>+ \ R J \E 


~d ) T f(E) 


(2) 


where, v is the frequency of the emitted photon, uj + and ui * 
are the statistical weights of the initial and final states re- 
spectively, f{E) is the free electron energy distribution func- 
tion, at is the photoionization cross-section for the inverse 
process of Equation[l]and R is the Rydberg energy constant. 
The other symbols have their usual meanings. In the present 
paper, the photoionization cross-sections for states of H and 


He are taken from Storey & Hummer ( 1991 1 and Hummer 


& Storey ( 1998 1 respectively as described in more detail by 


Ercolano & Storey ( 2006 ) . 

In this paper we consider two possible forms for f(E), 
the Maxwell-Boltzmann distribution 




(fcT) 


3/2 


E E 
— e kT , 
7 r 


(3) 


where T is the electron temperature, and the ^-distribution 
(|Vasyliunas||1968| |Summers fe Thorne||199l| given by 


A (e,t k ) = 


2 VE 


r (k + 1 ) 


§)lr( K - i) V ("-§ 


1 + ■ 


-(« + !) 


(4) 


where c is a parameter defining the distribution, Y is the 
gamma-function for the given arguments, and T K is a char- 
acteristic temperature. Note that f K becomes a MB distri- 
bution with T k — >■ T as k — \ oo. 

The normalization of the continuum flux to Hll flux re- 
quires that the effective recombination coefficients for Hll 
should be calculated both with a Maxwell-Boltzmann distri- 
bution and a ^-distribution. The hydrogen line emissivities 
tabulated by|Hummer &; Storey| |T987[) and|Storey fe Hum- 


1 19911 were calculated assuming that the free electron 


energy distribution is described by a Maxwell-Boltzmann 
distribution for all physical processes between bound and 
continuum states involving free electrons. Here we use the 
techniques and computer codes described in the last two 
references and extend them to include a K-distribution. At 
the electron number densities typical of photoionized nebu- 
lae (10 2 — 10 5 cm -3 ) the populations of the low-lying states 
of H are determined primarily by recombination and radia- 
tive cascading and are relatively insensitive to the ambient 
electron density. Hence any error introduced into the cal- 
culation of collision rates between high-n states caused by 
the use of a Maxwell-Boltzmann distribution rather than a 


^-distribution should have minimal effect at nebular densi- 
ties. We therefore make the approximation of computing the 
direct recombination coefficients to all the individual levels 
using a K-distribution but retain a Maxwell-Boltzmann dis- 
tribution for the energy and angular momentum changing 
collisions among the higher-n states. This should provide a 
good approximation for the Hll emissivity. 

For a general free electron energy distribution, the re- 
combination coefficient to a state X* is given by 


a RC — 


Rz ui’ 
y/2 c 2 m£ w 




On solving the collisional-radiative recombination problem 
for hydrogen we obtain the recombination coefficients to all 
levels and the effective recombination coefficients a e ff(A) for 
a transition of wavelength A. We define a line emission coef- 
ficient e(A) as the energy emitted per unit volume per unit 
time for unit ion and electron density, so that 

e(A) =a eff (A) (6) 


Nicholls et al (2013) derive the following approximate 


expression for converting recombination coefficients calcu- 
lated with a Maxwell-Boltzmann electron energy distribu- 
tion to those applicable with a ^-distribution, 


x(k)= Qk(A) 


£k(A) 


(!-^) r (« + l) 


«mb(A) £mb(A) (K-f)ir (k-A) 


(7) 


To obtain this result, |Nicholls et al (2013) use the fact 
that the photoionization cross-section falls approximately 
as (vo/is) 3 for v > t'o, where vq is the threshold frequency 
for photoionization, so that the integrand in Equation[5]con- 
tains a,l/v term. |Mcholls et a.Z[ (j2013 ) simplify the integral 
by moving the 1/v term outside of the integral, making it 
possible to carry out the integration analytically. In prac- 
tice the integral is a convolution of the free-electron energy 
distribution with the 1/v term, which is energy dependent. 
At low temperatures the rapid fall in the electron energy 
distribution function with increasing energy means that it 
is a good approximation to neglect the energy dependence 
of the frequency term. At higher temperatures, both terms 
are essential in the calculation of the recombination coeffi- 
cient. The expression for x{k) in Equation [7] is independent 
of temperature, density and transition. This approximation 
will fail at sufficiently high temperatures where the effects 
of the free-electron energy distribution and the frequency 
dependent term become comparable. 

Using the approximate function x(k), we can express 
the results of our more complete collisional-radiative treat- 
ment for e(A) in terms of a correction factor y( A, T, k) as 
follows 


i(A) = cmb (A) x(k) y{ A, T , «). 


(8) 


The values of £mb can be obtained from |Hummer fc Storeyl 
( 1987 1 and Storey & Hummer ( 1991 1. In Table[l]we tabulate 
values of y(X, T, k) for Hll corresponding to various values of 
data range of temperatures. We also tabulate in the last line 
i as a function of k. Note that, in principle, y also depends 
on electron density but in practice it is very insensitive to 
electron density for densities up to 10 5 cnY 3 , so we only 
tabulate values calculated at 10 3 cm' 3 . 
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Table 1. Values of the correction factor y(\,T,n), defined in 
Equation [8] as a function of re and logT [K] for computing the 
emission coefficient, e K (A) for Hll. In the last row of the table x. 
as defined in Equation [7] is given as a function of re. 


Table 2. Values of the correction factor z(2s,T, re), defined in 
Equation |10| as a function of re and logT[K] for computing the 
total recombination coefficient to the 2s state of hydrogen with 
re-distributed electron energies. 


logT [K] re logT [K] re 



50.0 

20.0 

10.0 

7.0 
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3.0 

2.0 


50.0 

20.0 

10.0 

7.0 

5.0 

3.0 

2.0 

2.6 

0.994 

0.994 

0.996 

0.997 

0.999 

1.004 

1.017 

2.6 

0.994 

0.994 

0.995 

0.995 

0.996 

1.000 

1.009 

2.7 

0.994 

0.995 

0.997 

0.999 

1.002 

1.011 

1.031 

2.7 

0.994 

0.994 

0.996 

0.997 

0.999 

1.005 

1.020 

2.8 

0.995 

0.996 

0.999 

1.001 

1.005 

1.017 

1.044 

2.8 

0.994 

0.995 

0.997 

0.999 

1.001 

1.009 

1.029 

2.9 

0.995 

0.997 

1.001 

1.004 

1.009 

1.024 

1.057 

2.9 

0.995 

0.996 

0.998 

1.000 

1.004 

1.014 

1.038 

3.0 

0.996 

0.998 

1.002 

1.007 

1.013 

1.031 

1.071 

3.0 

0.995 

0.996 

0.999 

1.002 

1.006 

1.018 

1.047 

3.1 

0.996 

0.999 

1.004 

1.009 

1.016 

1.038 

1.084 

3.1 

0.995 

0.996 

1.000 

1.003 

1.007 

1.022 

1.054 

3.2 

0.996 

1.000 

1.005 

1.011 

1.020 

1.045 

1.099 

3.2 

0.994 

0.996 

1.000 

1.003 

1.009 

1.025 

1.062 

3.3 

0.996 

1.000 

1.007 

1.013 

1.023 

1.053 

1.115 

3.3 

0.994 

0.996 

1.001 

1.004 

1.011 

1.029 

1.070 

3.4 

0.996 

1.000 

1.008 

1.016 

1.027 

1.062 

1.133 

3.4 

0.994 

0.996 

1.001 

1.005 

1.013 

1.033 

1.079 

3.5 

0.995 

1.001 

1.010 

1.019 

1.032 

1.072 

1.154 

3.5 

0.994 

0.997 

1.002 

1.007 

1.015 

1.038 

1.089 

3.6 

0.995 

1.001 

1.012 

1.022 

1.037 

1.083 

1.177 

3.6 

0.994 

0.998 

1.004 

1.009 

1.018 

1.045 

1.101 

3.7 

0.995 

1.001 

1.014 

1.026 

1.043 

1.096 

1.203 

3.7 

0.995 

0.999 

1.006 

1.012 

1.022 

1.052 

1.114 

3.8 

0.995 

1.003 

1.017 

1.031 

1.051 

1.112 

1.235 

3.8 

0.998 

1.001 

1.010 

1.017 

1.028 

1.062 

1.131 

3.9 

0.995 

1.004 

1.020 

1.036 

1.059 

1.128 

1.269 

3.9 

0.997 

1.002 

1.011 

1.019 

1.032 

1.069 

1.147 

4.0 

0.994 

1.004 

1.022 

1.040 

1.066 

1.144 

1.306 

4.0 

0.994 

1.000 

1.010 

1.019 

1.033 

1.075 

1.162 

4.1 

0.995 

1.006 

1.027 

1.047 

1.076 

1.165 

1.350 

4.1 

0.994 

1.000 

1.011 

1.022 

1.038 

1.085 

1.182 

4.2 

0.995 

1.008 

1.032 

1.054 

1.087 

1.187 

1.398 

4.2 

0.995 

1.002 

1.014 

1.026 

1.044 

1.096 

1.207 

4.3 

0.995 

1.010 

1.036 

1.060 

1.097 

1.210 

1.450 

4.3 

0.994 

1.001 

1.015 

1.028 

1.048 

1.107 

1.230 

x(k) 

1.008 

1.020 

1.043 

1.066 

1.103 

1.228 

1.596 










As can be seen, the values from Equation [7] are a good 
approximation for large values of re but the ratio increases 
as re decreases. In addition, the values from Equation [7] are 
generally smaller than the more exact values, with the dif- 
ference being largest for the smallest tabulated re and the 
highest temperatures, as expected. 

Since our model of the recombination process with re- 
distributed electron energies does not treat energy and an- 
gular momentum changing collisions correctly we cannot use 
the results to model the high Balmer lines as was done for 
example by Zhang et al ( 2004 1 using the code originally au- 
thored by one of us (PJS). We therefore restrict the current 
model to calculation of the continuum processes only, nor- 
malized to the Hll flux. 

Our model of the continuum also includes the hydrogen 
2-photon emission. The emission coefficient for the hydrogen 
2s-ls 2-photon emission is 


MB distribution in the same way as in Equation [8j via the 
function x(k) 

a«(2s) = q M b(2s) x(k ) z(2s,T, re) (10) 

where we tabulate the correction factors z(2s, T, re) in Ta- 
ble g] 

4 MODELING THE CONTINUUM 

The observed spectrum of Hf 2-2 with intensity normalized 
relative to Hll is shown in Figure [l] The figure also shows 
the spectral segments that we use for fitting, chosen to be 
as free as possible from significant spectral lines. The wave- 
length ranges for the three segments are A = 3585-3645, 
3780-3790 and 4180-4230 A. 

4.1 The continuum longward of the Balmer edge 


£2 q{v) = a(2s) 


A 2q {v) 


A 2q + Cb: 


s,2p 


Nr. 


hi/, 


(9) 


where a(2s) is the total recombination coefficient to the 2s 
state of H, A 2q is the total 2s-ls spontaneous transition prob- 
ability, A 2q {v) is the probability per unit frequency, C 2 a, 2 P 
is the coefficient for proton collisional transitions between 
the 2s and 2p states and N p is the proton number density. 
We assume that the population of the 2p state is negligible 
at the relevant densities (Hum mer fe Storey||1987| , and we 
take C' 2a , 2 p from Seaton (19551 and A 2q and A 2q (v) from 
Nussbaumer & Schmutz (1984 1 . Values of a(2s) for an MB 
distribution were provided by Hummer & Storey ( 1987 ) and 


Storey & Hummer ( 1991 1. From our new calculation of the 


hydrogen collisional-radiative problem with a re-distribution 
we obtain values of a (2s) as a function of re and tempera- 
ture. These values can be related to those obtained with an 


The observed continuum longward of the Balmer edge is 
comprised principally of hydrogen free-bound emission with 
principal quantum number n > 2, helium free-bound emis- 
sion, hydrogen 2-photon emission and scattered starlight 
from the central star of the nebula. The emissivity of the 
H 2-photon emission relative to Hll is relatively insensitive 
to the temperature of the emitting material and comprises 
up to 10% of the continuum in this spectral region. The con- 
tribution from n > 2 H free-bound emission varies strongly 
with the temperature, being 3% just longward of the Balmer 
edge at T = 10 4 K and falling rapidly as the temperature 
decreases. We model the remaining background contribu- 
tion to the flux relative to Hll with an empirical power-law 
distribution 


F( A) = F 0 



(ID 


© 0000 RAS, MNRAS 000 , 000-000 


The Continuum Emission Spectrum of Hf 2-2 5 


Table 3. Underlying continuum fitting parameters Fa and /3 
for two scenarios (a) A single MB distribution with tempera- 
ture 8800 K, and (b) a single re-distribution with re = 2 and 
T = 3000 K. 


FolA- 1 ] 

P 

a 0.161 

1.588 

b 0.174 

1.536 


where Fo is the contribution at the Balmer edge and wave- 
lengths in A are vacuum wavelengths. 

The parameters Fo and /3 were obtained by fitting the 
model continuum to the two longer wavelength segments 
shown in Figure[l]by minimizing the rms deviations of the fit 
from the observed spectrum. In computi ng the model co ntin- 
uum we take iV(He + )/./V(H + ) = 0.103 ( Liu et al|20 06). The 
same authors quote N (He 2+ ) / iV(H + ) = 0.002, which means 
that He 2+ recombination makes a negligible contribution to 
the continuum, so we neglect this component. The derived 
continuum is very weakly sensitive to the electron number 
density through the 2-photon component (Equation [9| and 
we adopt N e = 1000 cm -3 . 

The best fit values of Fo and /3 were derived for two 
models. The first uses a MB distribution for the free elec- 
trons with a temperature of 8800 K chosen to reflect the 


forbidden line values derived by Liu et al (2006). In any 


model comprising two MB distributions, the low tempera- 
ture component will contribute a negligible amount to the 
continuum longward of the Balmer edge, so it is appropri- 
ate to use the forbidden line temperature to determine the 
continuum component in this wavelength range. 

As we shall see below, the best fit to the continuum 
using a single re-distribution yields an electron temperature 
of a few thousand Kelvin and a small value of re. So for our 
second model of the underlying continuum we adopt a re- 
distribution with re = 2 and T = 3000 K. The results of 
the two underlying continuum models are summarized in 
Table [3] The values of Fo and /3 are relatively insensitive 
to the assumptions about the electron energy distribution 
and its temperature. To illustrate the fit longward of the 
Balmer edge, we show in Figure [l] the full fit to the two 
longer wavelength segments of the continuum including all 
the processes described above, using the parameters (a) from 
Table [3] and assuming a single MB distribution with T = 
8800 K. 


4.2 Fits with Maxwell-Boltzmann distributions 

It is clear from Figure [l] that the magnitude of the Balmer 
discontinuity and the slope of the Balmer continuum short- 
ward of the discontinuity cannot be explained by emission 
from hydrogen and helium at the forbidden line tempera- 


ture. This has already been discussed by Liu et al (2006) 


who estimated that the shape of the continuum implies that 
the emitting region has a temperature of « 900 K. In Fig- 
ure [2] we show the result of making a fit to all three spec- 
tral segments shown in Figure [I] using a single MB electron 
energy distribution in which the temperature is a free pa- 
rameter. The background continuum is computed using the 



Wavelength [Angstrom] 


Figure 2. The observed Hf 2-2 spectrum (solid) and the single- 
component Maxwell-Boltzmann fit (dashed) where the tempera- 
ture is treated as a free parameter. The final optimized tempera- 
ture is T = 1334 K. 


parameters from fit (a) in Table [ 3 ] The temperature of best 
fit is 1334 K and although the continuum model is improved 
compared to Figure [f] it matches neither the magnitude of 
the discontinuity nor the slope adequately. 

As discussed earlier, one resolution of the ORL/CEL 
abundance and temperature discrepancy in nebulae that 
has been proposed is that there are relatively cold high- 
metallicity knots embedded in hotter lower abundance gas. 
We therefore attempt to match the continuum with two 
components having different temperatures. We set the tem- 
perature of one component equal to a typical forbidden line 
temperature for this object of 8800 K and allow the other to 
vary. The continuum emissivity has only weak density de- 
pendence so we choose representative electron densities for 
the two components of 1000 cm -3 for the forbidden line re- 
gion and 5000 cm -3 for the cold material 
The relative emission measure of the two components is then 
allowed to vary by making the volume ratio of the two com- 
ponents a free parameter. In Figure |3]we show the resulting 
fit, which is excellent. The optimum temperature of the cold 
component is 540 K and the fraction of the total volume oc- 
cupied by the cold component is 0.00706. Hence the cold 
component has an emission measure, which is proportional 
to the product of the volume fraction and particle density 
squared, that is 15.0% of the total recombination line and 
continuum emission measure. 


(Liu et al 20061. 


4.3 Fit with re-distribution 

An alternative possible resolution of the ORL/CEL prob- 
lem is that there is no separation of the emitting material 
into high and low temperature components with significantly 
different abundances but rather a single medium in which 
the electron energy distribution is non-Maxwellian, with a 
re-distribution being a possible candidate. The Balmer con- 
tinuum close to the Balmer discontinuity provides a means 
of testing the validity of this proposal for the very low en- 
ergy part of the distribution function. In Figures [4] and [5] 
we illustrate how the magnitude of the Balmer disconti- 
nuity and shape of the Balmer continuum change as re is 
varied for two different temperatures. Compared to a MB, 
re-distributions have more particles at the lowest and highest 
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Wavelength [Angstrom] 

Figure 3. The observed Hf 2-2 spectrum (solid) and a two- 
component Maxwell-Boltzmann fit (dashed) with two tempera- 
tures: one is the forbidden line temperature which is fixed at 
T = 8800 K and the other is free to vary. The ratio of the vol- 
umes of the two components is also treated as a free parameter. 



Wavelength [Angstrom] 


Figure 4. Model continua with a K-distribution of electron en- 
ergies for a temperature T = 10000 K; k = 2 (dashed blue line), 
K = 5 (dotted green line), and ft = 25 (dot-dashed red line). 


energies and less at intermediate energies. At low energies 
the K-distribution function thus falls more rapidly with in- 
creasing energy than a MB distribution. This can be seen 
in Figure [5] for example, with the magnitude of the Balmer 
discontinuity increasing and the slope of the Balmer contin- 
uum becoming steeper as k decreases. 

We now fit the chosen three spectral segments to a k- 
distribution of which the temperature and value of k are 
free parameters and the final values that we obtained for 
these parameters are k = 2.11 and T = 4640 K. Figure [6] 
shows a plot of this fit. On visual inspection the best-fit k 
distribution where k and temperature are allowed to vary 
is significantly less good than the two-component MB fit. 
It should be remarked that the optimized k obtained from 
this fitting exercise is very different to values proposed by 


Nicholls et al (2012) which usually fall in the range 10-20. 


For comparison, we include in Figure[6]a fit carried out with 
a fixed value of ft = 10. The temperature was allowed to vary 
and the best fit value was 1665 K. In the next section we 
attempt to quantify the relative quality of these fits. 



Wavelength [Angstrom] 


Figure 5. Model continua with a fs-distribution of electron en- 
ergies for a temperature T = 1000 K; ft = 2 (dashed blue line), 
ft = 5 (dotted green line), and ft = 25 (dot-dashed red line). 



Wavelength [Angstrom] 


Figure 6. Fits obtained by simultaneous optimization of ft and 
T (dashed red line) and by fixing ft = 10 and allowing T to vary 
(dotted blue line). 


4.4 Statistical Analysis 

We wish to address the question of whether the fit to the con- 
tinuum using a fc-distribution is significantly worse than the 
two-component MB fit. A visual inspection of Figures [3] and 
[6] shows that while the-two-component fit provides a very 
good match to the three segments of the continuum used for 
fitting, the best ^-distribution fit underestimates the mag- 
nitude of the Balmer discontinuity and also underestimates 
the steepness of the decline at shorter wavelengths. We will 
use a comparison of \ 2 values for these two fits to quantify 
the relative goodness of fit. To estimate \ 2 we need an esti- 
mate of the uncertainties on the observational data, which 
by inspection and from instrumental considerations are not 
independent of wavelength, being much larger at the short- 
est wavelengths. We therefore compute separate rms devia- 
tions of the data from the two-component MB fit for each 
of the three wavelength segments. Since we use deviations 
from the two-component fit to define the data uncertainties 
it follows that the value of y 2 for the fit, summed over the 
three wavelength segments, will be given by 

X 2 = N df = N dp - N fp = 237, (12) 
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where Ndf is the number of degrees of freedom, Nd p is the 
number of data points in the three segments and Nf p is 
the number of free parameters. We now use the estimates 
of the data uncertainties from the two-component MB fit 
to compute \ 2 f° r the best fit K-distribution and we find 
X 2 = 433. We wish to test the null hypothesis that the k- 
distribution is an equally good fit to the data as the two- 
component MB fit. The fol l owing squa re-root transforma- 
tion of \ 2 (Abramowitz fe Stegunji972 1, 


t=^2^~ sj 2Ndf , 


(13) 


is such that for large Ndf, the transformed distribution is 
approximately normal with zero mean and unit standard 
deviation. For the best fit K-distribution we find that t = 
7.67 standard deviations from the zero mean, indicating that 
the probability that the null hypothesis is true is vanishingly 
small. 

Thus, a two-component Maxwell-Boltzmann model is 
highly favored compared to the single-component k model. 
The single-component MB model is even less likely to be 
correct than either of the two-parameter models based on 
the same \ 2 analysis, giving t = 14.4 standard deviations. 


5 CONCLUSIONS 

The analysis of the Balmer continuum spectrum of Hf 2-2 
seems to indicate the signature of a two-component nebula 
with two different temperatures. Assuming that one compo- 
nent has a temperature of 8800 K, we find that the best fit 
to the data is obtained by adding a second component with 
temperature of 540 K and emission measure 15.0% of the 
total. We also modeled the continuum with a single com- 
ponent with K-distributed free electron energies, but on the 
basis of a \ 2 test we found this model to be significantly 
less likely to be correct. In the case of K-distribution mod- 
els, it is important to note that modeling the Balmer con- 
tinuum only samples the free electron energy distribution at 
the lowest energies, below the mean energy, and therefore 
gives no information about departures from MB at energies 
above the mean energy which might affect the excitation of 
nebular forbidden lines, for example. The extreme nature 
of Hf 2-2 and the similarity of some of its spectral features 
to old novae like DQ Her may also cast some doubt on the 
applicability of the conclusions reached to the properties of 
planetary nebulae in general. 
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